Off equilibrium dynamics of the Frustrated Ising Lattice Gas 
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We study by means of Monte Carlo simulations the off equilibrium properties of a model glass, 
the Frustrated Ising Lattice Gas (FILG) in three dimensions. We have computed typical two times 
quantities, like density-density autocorrelations and the autocorrelation of internal degrees of free- 
dom. We find an aging scenario particularly interesting in the case of the density autocorrelations 
in real space which is very reminiscent of spin glass phenomenology. While this model captures the 
essential features of structural glass dynamics, its analogy with spin glasses may bring the possibility 
of its complete description using the tools developed in spin glass theory. 



Much effort is being currently devoted to reach a 
reasonable theoretical understanding of structural glass 
physics. While the huge amount of experimental work 
available shows an extremely rich phenomenology, a suc- 
cesfuU description in terms of microscopic models is still 
lacking. The theoretical models available are mainly phe- 
nomenological, basically focusing the non-Arrhenius re- 
laxation in the so called fragile glasses jl). Among the 
most succesfull ones are the free volume model || and 
the entropy model || in which the relevant variables for 
the description of a glass transition from a supercooled 
liquid phase are respectively volume and entropy. While 
the last model predicts the existence of a thermodynamic 
second order transition at a well defined temperature, it 
is extremely difficult to obtain evidences of it as it is in 
practice experimentally inaccessible. Another succesfull 
approach has been the Mode Coupling Theory of Gotze 
and Sjogren ph. This dynamical approach is qualita- 
tively correct in predicting the behavior of time corre- 
lations and responses in supercooled liquids. Recently 
is has been shown to correspond to the high temper- 
ature limit of a dynamical theory of spin glasses ||. 
This raised the interesting possibility of describing the 
structural glass physics by exploiting the analogy with 
some spin glasses in which the transition is discontinuous 
as originally suggested by Kirkpatrick et al. Q. More- 
over, some finite dimensional models have been proposed 
whose mean field limit would provide such a transition 
JjJ. Nevertheless, although spin and structural glasses 
have some basic features in common, like the character- 
istic slow dynamics, an essential difference is the absence 
(presence) of quenched disorder in structural glasses (spin 
glasses). Unlike in spin glasses, the nature of the glass 
transition may be purely dynamical in origin, without an 
underlying thermodynamic transition. The existence of 
a growing correlation length, for example, has not been 
considered up to very recently in structural glass mod- 
els H . The absence of a simple microscopic model turns 
these questions very difficult to be answered and much of 
the present knowledge comes from computer simulations 



of Lennard-Jones systems S or from simple kinetically 
constrained models jlO] ^3|- These, in particular, have 
dynamically self-induced frustration effects obtained by 
restricting the possible Monte Carlo movements and re- 
produce reasonably well the glassy phenomenology. How- 
ever, although having the advantage of being lattice mod- 
els, it is not obvious how to relate the ad hoc kinetic rules 
with the underlying physics. 

A possible candidate to fill in this gap is the Frus- 
trated Ising Lattice Gas model (FILG) (l4| , a Hamilto- 
nian lattice gas in which the presence of internal degrees 
of freedom subjected to quenched disorder mimics the 
geometric frustration that slows down the motion of the 
molecules as the system is cooled or compressed. The 
introduction of internal degrees of freedom is responsible 
for the slowing down of the diffusional dynamics of the 
particles and the necessity of considering them in order 
to describe the essence of glass physics has recently been 
addressed by Tanaka p5[ . The thermodynamic proper- 
ties of the FILG in 3D as well as its equilibrium dy- 
namics have been studied by Nicodemi and Coniglio [|l6| . 
The model shows many glass properties, including depen- 
dence on the cooling (or compression) rate, stretched ex- 
ponential behavior in correlation functions, a dynamical 
singularity in which the diffusion constant goes to zero, 
the breakdown of the Stokes-Einstein relation along with 
anomalous diffusion at intermediate times. As in many 
glasses no singularities in the linear susceptibilities are 
observed, in particular the compressibility is continuous 
everywhere. Scarpetta et al. pi studied an equivalent, 
non local version of the FILG (the Site Frustrated Per- 
colation, SFP) in 2D, finding a behavior similar to the 
3D FILG, the main difference being the (Arrhenius) dy- 
namical singularity occurring at zero temperature. This 
model also seems to relate the glass transition with a 
percolation-type transition as evidenced by the onset of 
several precursor phenomena |l7j . 

Non equilibrium phenomena, aging being one example, 
are widespread on a great variety of systems, including 
polymers, granular materials and ferromagnetic coarsen- 
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ing ||l9|-^H and appears in the relaxation dynamics of 
both spin and structural glasses, either theoretically and 
experimentally. In this letter we report results for the off 
equilibrium dynamics of the FILG in three dimensions, 
showing a characteristic aging dynamics present in some 
two times quantities and some possible scaling scenarios. 
Of particular interest are the results for the density auto- 
correlations which suggest the definition of a non-linear 
compressibility which may bring information on a pos- 
sible thermodynamic transition analog to the magnetic 
transition in spin glasses. 

The FILG is defined by the Hamiltonian: 

H = — J {sijcricrj - l)mnj - n,-. (1) 

<ij> i 

There are two kinds of dynamical variables: the site oc- 
cupation m = 0, 1 (i — 1 . . . N) and the particles internal 
degrees of freedom, oi = ±1. The usually complex spa- 
tial structure of the molecules of glass forming liquids, 
which can assume several spatial orientations, is in part 
responsible for the geometric constraints on their mo- 
bility. Here we take the simplest case of two possible 
orientations, and the steric effects imposed on a particle 
by its neighbors are felt as restrictions on its orientation 
due to the quenched random variables = ±1. The 
key role of the first term of the Hamiltonian is that when 
J — > oo (recovering the SFP) no frustrated link can be 
fully occupied, implying that any frustrated loop in the 
lattice will have a hole and then p < 1 preventing the 
system from reaching the close packed configuration. Fi- 
nally p represents a chemical potential ruling the system 
density (at fixed volume) and, by taking p — > oo we re- 
cover the Edwards- Anderson spin glass model. It should 
also be mentioned that, after including gravity, the same 
model successfully describes granular materials under vi- 
bration fH,|23| > another class of systems where geometric 
frustration rules its behavior. 

By increasing p the model presents two characteristic 
points |l^,[l7]]. For p < 0.75 (low density) it shows liq- 
uid like behavior, time correlation functions decay expo- 
nentially, equilibration is quickly achieved and the parti- 
cles mean squared displacement grows linearly with time, 
a simple diffusion scenario. At p « 0.75 there is a 
percolation transition (the corresponding density being 
p Ri 0.38). Dynamically it manifests in the onset of two 
different relaxation regimes in the correlation functions, 
a fast exponential relaxation at short times and a slow 
relaxation at longer times characterized by stretched ex- 
ponentials. The diffusion is still linear for long times but 
the diffusion coefficient becomes smaller as the density 
grows. Also, for a fixed p, the equilibrium density de- 
pends on the cooling rate. The dynamics becomes slower 
as the chemical potential grows (or equivalently, the tem- 
perature is lowered) and a second transition is reached for 
p > 6. This is a spin glass transition associated with the 
frozen in of the internal degrees of freedom. Interestingly, 
at this point a dynamical singularity is also present, man- 
ifested by the vanishing of the diffusion constant. Besides 



the qualitative analogies with the physical processes typ- 
ical of glass forming liquids, there are a few points of 
more fundamental character. The presence of two char- 
acteristic temperatures (or chemical potentials) separat- 
ing different dynamical regimes is common with a class 
of mean field theories of spin glasses and also with what 
is observed in real glass formers. The percolation tran- 
sition corresponds to the dynamical transition in mean 
field p-spin or Potts glasses or with the glass transition 
in the mode coupling theory. In the FILG the dynamical 
signature of the transition is consequence of the appear- 
ance of percolating clusters, a geometrical feature which 
may be responsible for the slowing down of the dynamics 
in structural glasses. The second transition corresponds 
to the ideal glass transition at which the relaxation times 
diverge and the diffusion constant goes to zero. In the 
FILG the structural manifestation of this transition is 
the presence of a frozen percolating cluster. The density 
of the system approaches a critical value p c « 0.7, at- 
tainable only by an infinitely slow cooling. These facts 
turn the FILG an important model for studying different 
aspects of the glass transition. 

Different protocols can be envisaged for studying, for 
example, correlations and response functions. We have 
prepared the model in a non equilibrium state, setting the 
parameters as in |]l6| with an initial density lower than 
the critical one. Then the system is quenched to a su- 
per critical chemical potential and dynamic correlations 
are recorded. We have used a Monte Carlo dynamics 
which alternates flipping of the internal degrees of free- 
dom, creation and destruction of particles in a plane sur- 
face (which mimics a compression experiment) as in |ll]] 
and particle diffusion. The creation-destruction of parti- 
cles is fast enough as to destroy all particles at long times 
if it is allowed in the bulk. 

Consider the connected two point correlation 

c(t, t w ) = — ni(t + t w )n 2 (t w ) - p(t + t w )p(t w ) (2) 

i 

where the global density at time t is given by p(t) = 
N~ 1 ^2 i rii(t). We now define the density autocorrela- 
tions as C n (t,t w ) — c(t,t w )/c(0,t w ). In figure (|]) the 
behavior of C n (t,t w ) is shown as a function of t in a 
semi-log plot after a quench in chemical potential to a 
value p = 10, for waiting times between 2 5 and 2 17 . 
A typical aging scenario is present signalling the slow- 
ing down of the dynamics as the waiting time grows. 
For the longest waiting times the correlation presents a 
rather fast relaxation to a plateau in which the system 
is in quasi-equilibrium: the dynamics is stationary and 
the fluctuation-dissipation relations hold. The plateau 
separates two time scales typical of glassy systems: a (3 
(fast) relaxation for small time and an a (slow) relax- 
ation at longer times, corresponding respectively to the 
fast movements of the particle inside the dynamical cages 
and the large scale, cooperative process that takes much 
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more time in order to rearrange the cages. Moreover, in 
this very long time regime (t 3> t w ), the system falls out 
of equilibrium, the correlations decay to zero asymptot- 
ically and time translational invariance (TTI) no longer 
holds with the corresponding violation of the fluctuation- 
dissipation theorem (FDT). 




LJ ' "J ' "J ' "J ' "J ' "J 1 

1 10 100 1000 10000 100000 

t 



FIG. 1. Density autocorrelations after a quench to n = 10 
at T = 1 and J = 10 for L = 20. The waiting times range 
from 2 5 (bottom) to 2 17 (top) and the averages are over 50 
samples. 
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FIG. 2. Scaling plot of the five upper curves of figure (|l|). 
The initial five points (/3 decay) were discarded and the pa- 
rameters are A — 0.59, a — 0.46 and /i = 0.9. See text for 
details. 

As the C n (t,t w ) curves have a complex behavior it is 
possible to scale the curves only after a careful analysis of 
the different time scales present 0] . It is clear that only 
for the largest waiting times different relaxation scales 
can be observed. We restrict our analysis to the five up- 
per curves (t w = 2 13 . . . 2 17 ). In these curves, a fast initial 
decay is observed up to times t ~ 100 MCS, correspond- 
ing to the already mentioned [3 decay which we will not 
analyze in any more detail. Then a plateau characteris- 
tic of quasi-equilibrium dynamics develops. Its asymp- 
totic value as t w — > oo defines the ergodicity breaking 
(or Edwards- Anderson) parameter. For a fixed, large t w 
and t 3> t w , the system begins to fall out of equilibrium 
and the correlations decay to zero. But only the very 
early epochs of this decay can be observed within the 



times of our simulations, although we can observe the 
crossover between equilibrium and non equilibrium dy- 
namics. Consequently, in order to obtain a good scaling 
for this last regime we have subtracted from C n (t,t w ) a 
stationary contribution of the form: 

C n (t,t w ) = Coo + At- a , t w &xed,t<t w . (3) 

with Coo , A and a as fit parameters. For the non sta- 
tionary regime we assumed a time dependence of the 
form h(t w )/h(t + t w ) with h(x) given by Q h(x) = 
exp [(1 — /i) _1 (x/r) 1_A1 ] with \i < 1 and r a microscopic 
timescale. Note that this form is quite general as one re- 
covers the cases of a simple t/t w dependence (full aging) 
when fi = 1 and stationary dynamics when /i = 0. The 
final scaling is shown in figure (||). The scaling obtained 
is much better than assuming only full aging or activated 
dynamics even after including in these cases the contri- 
bution from the stationary decay. 

We have also measured, as shown in figure the 
autocorrelations of the internal degrees of freedom: 

C s (t,t w ) = J^^2 s i{t + t w )n,i(t + t w )si{t w )rii{t w ), (4) 

i 

The internal degrees of freedom correspond to a diluted 
Edwards- Anderson (EA) spin glass and the curves should 
be compared to the ones of the 3D EA model |2(| . From 
our knowledge of the 3D EA model one may expect the 
scaling of the autocorrelation in the aging regime be 
of the form t~ a C(t/t w ). Nevertheless we verified that 
an activated dynamics scaling of the form C s (t,t w ) ex 
C (\og(t) / \og(t w )) works better. Moreover, doing an 
analysis similar to the one for C n (t,t w ) the scaling can 
be slightly improved (see fig.(||), inset). 
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FIG. 3. Autocorrelation of the internal degrees of freedom 
after a quench to fi = 10 with T = 1 and J = 10 for L — 20. 
Waiting times from 2 s to 2 17 (top to bottom) and average over 
50 samples. Inset: scaling with to = 0.5, /i = 0.9, a — 0.091 
and A = 0.057 (see text). 

To our knowledge this is the first Hamiltonian lattice 
model which presents the essentials of structural glass 
phenomenology. In particular the results for the density 
autocorrelation are very promising suggesting that it may 
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be possible to apply spin glass ideas and techniques for 
an analytical investigation of the model |L8|]. The simi- 
larity between the density autocorrelation and the spin 
glass overlap function suggests the introduction of a non- 
linear compressibility in analogy with the spin glass sus- 
ceptibility, which may be a good quantity for studying 
the possibility of an underlying phase transition with a 
growing correlation length in the model. Another issue 
that must be studied is the precise form of the violation 
of FDT through the so called fluctuation dissipation ra- 
tio p7| . Also other protocols may be implemented in 
order to probe the off equilibrium dynamics of the sys- 
tem. Work in progress, keeping fixed the global density 
and performing a quench in temperature, indicates that 
the phenomenology is similar to what we have presented 
here. Moreover, the evaluation of the root mean square 
deviation of the particles and the incoherent scattering 
function (which is related to the Fourier transform of the 
density correlations) also show evidence of slow dynam- 
ics and aging. The results of these investigations will be 
published in a future paper. We are working also in the 
2D version of the model: here, some differences are ex- 
pected with respect to 3-D because the dynamical singu- 
larity only occurs for T = (/i — > oo) and the relaxation 
time diverges with an Arrhenius law. Still other issues 
can also be explored. For instance, to what extent, if any, 
the scenario for the 3D model presented here changes as 
one goes to infinite range connections. In the mean field 
version, where both first and second order transitions 
show up, it might be possible that different aging regimes 
are present Jig ] , and would be interesting to know if the 
FILG is a finite dimensional version of a model whose 
mean field limit has a discontinuous transition. 
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